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Abstract 

We revisit the classical approach of comoving coordinates in relativistic hy- 
drodynamics and we give a constructive proof for their global existence under 
suitable conditions which is proper for stochastic quantization. We show that 
it is possible to assign stochastic kinematics for the free relativistic spinless 
particle as a Markov diffusion globally defined on M 4 . Then introducing dy- 
namics by means of a stochastic variational principle with Einstein's action, 
we are lead to positive-energy solutions of Klein-Gordon equation. The pro- 
cedure exhibits relativistic covariance properties. 
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I. INTRODUCTION 



Probability Theory, and Markov processes in particular, are currently recognized as key 
elements for a proper modelization of various relativistic phenomena within the framework 
of the Euclidean formulation of Quantum Field Theory. On the other hand, the construction 
of a mathematically consistent description of diffusive phenomena in Minkowskian spaces 
must face severe problems. This is due essentially to the difficulty to make the Markovian 
property compatible with the required relativistic covariance (for different approaches and 
proposals concerning this longstanding open problem see for example Refs. |l|-|8j and ||). 
A satisfactory theory seems to be still lacking. 

The aim of this work is to introduce a new method for treating phenomena described 
by Markovian diffusions in the configurational space within a relativistic setting. As a first 
application we procede to perform a covariant stochastic quantization of the free relativistic 
spinless particle. From a physical point of view, the idea underlying our construction can be 
sketched as follows. Let us consider, for the sake of simplicity, a three-dimensional Markov 
diffusion with constant coefficient equal to u, satisfying Ito's stochastic differential equation 

dq(t) = b(q(t),t)dt + v* dw(t) , te[0,+oo), (1.1) 

where b is a drift-field and w a standard Wiener process. Introducing the time dependent 
density p of the process q(t), one can show that under some regularity assumption |T0|JTT[| 
there exists a velocity field v (usually called the "current velocity field") such that the 
following continuity equation holds: 

d t p + V-(pv) = 0. (1.2) 

Thus, if b is sufficiently smooth, to every Markov process given by a solution of ( |1.1| ) can 
be associated an "hydrodynamical (Eulerian) structure" represented by the time dependent 
couple (p,v), which satisfy the constraint (|1.2|). 

Let us suppose that, together with the diffusion, we are interested to consider, at a 
generic point q and time t, basic observable quantities as energy-momentum tensor, entropy, 
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temperature and so on. Thus, in agreement with the approach usually adopted in the 
relativistic hydrodynamics of a perfect fluid, we would define those objects in the frame of 
an observer which moves at point q and time t with velocity v(q, t) (the "comoving observer" 
at point q and time t). Intuitively, it is therefore reasonable to extend this prescription also 
to the "probabilistic" objects we are considering. This idea was in fact already exploited in 
Ref. fl2| to show that it is possible to redefine Nelson's stochastic kinematics in a "small" 



neighborhood of properly defined "comoving observers"; starting from this and a suitable 
extension of the Einstein action, the Klein-Gordon equation for the free relativistic particle 
was there derived without ad hoc assumptions. 

In the present work we will go beyond the local approach outlined above, and we in- 
troduce instead a "global comoving coordinate system". A non linear transformation is 
performed on flat space from natural coordinates to a class of comoving ones. Our discus- 
sion will not depend on the particular choice in this class and the results will be expressed 
in the language of general curvilinear coordinates. The method is intimately connected with 
classical results as Frobenius theorem for vector fields and it is close to the construction of 
cosmic standard coordinates in General Relativity (see for example Ref. [HJ). We will give 
a precise and simple formulation of the basic mathematical facts which does not seem to be 
commonly available in the literature and which, in particular, is proper for the application 
to stochastic quantization. 

The content of the work is organized as follows: in Sec. II we give a sufficient condition for 
a vector field V on R m , m > 1, induces a "global comoving coordinate system". In Sec. Ill 
the same result is reformulated for vector fields on M 4 and pseudo-Riemannian manifolds. 
In Sec. IV we apply these results to redefine Nelson's stochastic kinematics for a free 
relativistic spinless particle and to derive in a simple and fully consistent mathematical way 
the Klein-Gordon equation. The stochastic interpretation of different subsets of solutions to 
this equation is given in Sec. V, where the non relativistic limit is also described. Finally, 
in Sec. VI we discuss some possible physical implications and briefly consider some future 
developments and applications. 



II. VECTOR FIELDS GENERATING GLOBAL COMOVING COORDINATES 

IN EUCLIDEAN SPACES 



Let us consider a vector field V : R 1+n — >■ R 1+ri , x t— > V(x). We shall use for the 
components of x and V the notation x = (x°, x 1 , . . . , x n ) and V = (V°, V 1 , . . . , V n ). We 
are searching conditions on V so that there exists a coordinate transformation <3> which 
changes the algebric vector V(x) into V'(x) := (|| V(a;)||, 0, . . . , 0) for every x G R 1+n . The 
inverse function theorem allows to prove the following: 

Theorem 1: Let V be a vector field on R 1+n with components (V , V 1 , . . . , V n ), which 
satisfies the following conditions: 

i) V°(x)^0 Vi 6 R 1+n , 

ii) There exists a scalar field S G C r (R 1+n , R) ,r > 2, such that 

V(x) = S(x) , V/i = 0, 1, . . . , n , Vx E R 1+n . 
Then there exists a diffeomorphism $ : R 1+n — > R 1+n such that 

JL d^(x) \\V(x)\\ if a = , 

^=o I otherwise . 

Proof: Since is C 1 , then there exists the congruence Cy generated by V. That is, given 
any point x G R 1+ri , there exists a unique integral curve j x G Cy to which x belongs. For a 
fixed xo G R 1+ ™ let us consider the equation 

S(x) - S(x ) = . (2.1) 

Denoting thus yo the n-dimensional vector (xq, . . . ,Xq), by assumptions i) and ii) and 
the implicit function theorem there exists a neighborhood I yo C R™ of and a function 
/ G C^/yo, R) such that 



S(f(x 1 ,...,x n ),x 1 ,...,x n ) - S(x )=0 y(x 1 ,...,x n )el 



yo ■ 



Let now T denote the border of I yo and p any point belonging to T. Then the lim q _ >p /(q) 
may be finite or infinite (the existence is ensured by the fact that the application of the 



implicit function theorem in p would otherwise lead to a contradiction). In the first case 
one can apply again the implicit function theorem at p, and so on. Thus for each xq G R 1+rt 
equation ( |2.1| ) defines a unique regular n-dimensional hypersurface S Xo in R 1+n . We also 
observe that the vector 



(dS(x)_ dS(x)_ dS(x) \ 
{ ' ' \ dx° ' dx 1 dx n J 



is by construction orthogonal to the surface in the point x. Hence we can conclude that for 
each point Xq in R 1+n there exists a unique hypersurface E^ and a unique integral curve 
7 Xo to which Xq belongs. In addition, 7^ is orthogonal to £ Xo in x - In other words, the 
congruence Cy is globally hypersurface-orthogonal. 

To construct the new coordinate system in R 1+n let us fix a point O' as new origin 
and consider 70' (with the natural orientation) and Eo'- Being Eo' a regular hypersurface 
imbedded in R 1+n , it is also a n-dimensional Riemannian manifold. Let us introduce on 
it a coordinate system q 1 ,. . . ,q n (notice that this surface has a unique chart R n ) and let 
(jjj , i,j = 1, . . . , n be its metric tensor, induced by the Euclidean scalar product in R 1+ ". 
Then to every point x G R 1+n there corresponds a unique point on 70' (intersection between 
70' and E x ) with coordinate X(x) and a unique point on So/ (intersection of Eo' and 7 2 ) 
with coordinates q ,1 (x), . . . , q n (x), and viceversa. Therefore the map (see Fig. 1) 

$ : R 1+,t — ► R 1+n , 

x ^ ^.= (X(x),q 1 (x),...,q n (x)) (2.2) 

is one-to-one. Existence and continuity of first-order partial derivatives of $ immediately 
follows from the regularity of the integral curves and of the hypersurfaces we are considering. 

Let us now denote (p x the solution of the ordinary differential equation 
with initial condition 



if x (r) = x . 

Hence we have 

de(r) = ^ 7 d V :(r) + 0(\\d Vx (r)\\ 2 ) 



dx 

= ^ 7 V»(x)dT + 0(dT 2 ). (2.3) 

Let us now choose X(x) equal to the length of the arc on j ' an d £° := A(x). Then we 
observe that defining V'^ by the equality 

V»==^V, (2.4) 

we immediately see that, since d^{r) = for /x = 1, . . . , n and d£°(r) = \\V(x)\\ dr + 0(dT 2 ), 
only the component V'° is different from zero and it is equal to ||V(a;)||. □ 

It is not necessary to particularize the coordinate A to be the length of the arc on 70'. 
In fact we can prove the following: 

Corollary 1: Let £° = A be any coordinate on 70'. Introducing on jo' the positive metric 
function /(A) we have 

-WL if m = o, 

otherwise . 

Proof: Denoting by l(X) the length of the arc (O', A), we can rephrase eqs. fl2.3Q - (|2.4fl by 
observing that 

<%<> = dX= = in^fl dr + ( dr 2) = V'°(0dT + O(dT 2 ) , 

from which the result follows. □ 



The coordinate transformation $ transforms the Euclidean metric tensor 5^ into a new 
one, which we denote by g^ v . This transformed tensor has in turn a particular ily simple 
structure, as one can immediately recognize by observing that any vector dx G R 1+ ™ applied 
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at the point x can be decomposed into two orthogonal vectors dz and dy, where dz is tangent 
to 7 S at x and dy is tangent to S s at x. Thus we can write 

\\dx\\ 2 = \\dz\\ 2 + \\dy\\ 2 

= f(\)d\(x) 2 + jr a ij (q 1 (x),...,q n (x))dq i (x)dq^x)+0(\\dx\\ 4 ), 

where as before /(A) and are the components of the metric tensor of 70' and S G / in the 
given coordinate \(x) and q 1 (x), . . . , g n (x) respectively. Thus we can state the following: 

Corollary 2: The global coordinate transformation $ : x £ := $(:r) , x G R 1+ ™, 
transforms the Euclidean metric tensor 5^ into the tensor g E v with components 



fj, — , V — 1, ... ,n ; V — , fj, — 1, ... ,71 , 

<V(£\ ■ ■ ■ ,t, n ) otherwise. 



Owing to the explicit structure of the transformed metric tensor, we are able to com- 
pletely characterize our construction from the geometric point of view. In particular, being 
the map $ a global coordinate transformation on flat space 

R l+n 

we know that the Riemann 

curvature tensor associated to g E identically vanishes. We shall see now that this property 
extends to a, so that £0' eventually has only an extrinsic curvature. We summarize the 
result in the following: 

Corollary 3: The hypersurface £0' has zero intrinsic curvature. 

Proof: We need to prove that the n-dimensional Riemann curvature tensor R^im , 
i,k,l,m — 1, ... ,n associated to a is equal to zero. Starting from the (1 + n)-dimensional 
tensor associated to g E , let us express it, as usual, in terms of Christoffel symbols: 

R Mm = ^Tp ~~ ~Q£ m ~ + ^al^tm ~~ ^am^tl a = 0, . . . , 1 + 71 . (2.5) 

From Corollary 1 and the well known general relation 



one can easily check that the connection T induced from g E has a particularity simple 
structure, which we write as 

r = r (1) © r (n) , 

with I^ 1 ) = r[} and = {Tj k }, i, j, k = 1, . . . ,n corresponding to the metrics on 
and S / respectively. Exploiting this fact, eq. (|2.5|) becomes 

R (l+n)i _ R (n)t , pjpO _ -pi r _ p(»)' 
rt klm ~ rL klm' 1 0l 1 km 1 Om 1 kl ~ " klm ■ 

The left hand side is zero by flatness, so that necessarily R^^ = 0. □ 



It is also worth to stress that, as usually done in the theory of Riemannian manifolds, we 
can introduce normal coordinates . . . , Q in a neighborhood I x of each point x G R 1+n . 

This is obtained by choosing normal coordinates both in the one- dimensional neighborhood 
of the projection of x on 70' and in the n-dimensional neighborhood of the projection of x on 
Eq/ along Cy. (Indeed, due to the flatness of S / (Corollary 3), this can be done globally, 
so that it possible to define a coordinate system that yields an orthonormal coordinate 
basis at every point of R 1+n ). Hence, in case O' coincides with O, the linear part of the 
transformation $* : y 1— > := $*(y) can be thought as a rotation on R 1+n which aligns the 
x°-axis along the direction of the vector field and which is different from point to point if 
the latter is not constant. 



III. GLOBAL COMOVING COORDINATES IN MINKOWSKIAN SPACES AND 

PSEUDO-RIEMANNIAN MANIFOLDS 

In this section the Minkowski spacetime M 4 is understood to be endowed with the stan- 
dard Euclidean topology. (By the way, in the literature the so-called "path topology" has 
also been considered ||14|| , nevertheless it looks quite unsuitable in our setting). Then it is 



straigthforward to extend Theorem 1: 



S 



Theorem 1': Let M 4 be equipped with the Euclidean topology and let rj^ := 
diag(— 1, 1, 1, 1) the Minkowskian metric tensor. If V is a vector field on M 4 with covariant 
components (Vq, Vi, V2, V3) satisfying the following conditions: 

i) V {x)^0 Vi G M 4 , 

ii) There exists a scalar field S G C r (M 4 , R), r > 2, such that 

V,(x) = ^S(x) , V/i = 0, . . . , 3 , Vx 6 M 4 , 
then there exists a diffeomorphism $ such that 



u=0 UX 



V^x)V»(x) I if n = , 

otherwise . 



Proof: With the given topology, M 4 is a pseudo-Riemannian manifold with a unique 
chart R 4 , cartesian natural coordinates and metric tensor 77^. Thus the proof of Theorem 

1 remains the same as far as it does not involve the metric tensor. This occurs only where 
we assert that the vector V(x) is orthogonal to T, x at point x. But we have seen that the 
Euclidean product between a vector with (covariant) components dS(x)/dx^ and any vector 
e(x) G R 4 which is tangent to S G / at point x is equal to zero. Hence we can write 

JL- S ( x )e»(x) = V^x)e»(x) = v^V^e^x) = , 

so that, by the arbitrariness of the tangent vector e(x), the vector V(x) = {V°, V 1 , V 2 , V 3 ) 
is orthogonal to T, x at the point x in the Minkowskian metric. □ 

We shall now for definiteness suppose that the congruence C v is timelike. Corollaries 1, 

2 and 3 can be easily extended in the following way. 

Corollary 1' : Let £° = A be any coordinate on (0',\). Then introducing on 70' the 
negative metric function g 00 (A) we have 

' y/-V»{x)Vp{x) 



V'^ix) 



, if ji — , 

V-Soo(A(rr)) 

otherwise . 



Corollary 2': The global coordinate transformation $ : x i— > £ := $(x) , £ G M 4 
transforms the Minkowskian metric tensor 77^ into the tensor g^ v defined by 





<v(e\e a ,e 8 ) 

Proof: This structure immediately comes from that of g E . In fact we can write 



\i = v = , 

fi = 0,i/=l,2,3;i/ = 0,fi = l,2,3 
otherwise . 



9%{Q = <W" 



-(0 (0 



and 



(3.1) 



(3.2) 



so that a comparison between ( |3.1| ) and ( |3.2| ) yields g^ v = g E for /x 7^ and f 7^ while 

rf) = = -/(a- □ 

Corollary 3 still holds (the same proof with n = 3). 



Exactly as before, it is possible to introduce normal coordinates in a neighborhood 
of any point x G M 4 by projecting 7^. on So' along Cy. In this case, if O' coincides with 
O, the coordinate transformation : y 1— > := $*(?/) , y E I x has the linear part which 
reduces, up to a spatial rotation, to a Lorentz boost associated to V(sc), or, explicitely: 



V da; 1. 



(x) = A-^Cx) 



Vi G M 4 . 



Theorem 1' can also be generalized to vector fields defined on a pseudo-Riemannian 
manifold M. In fact, given a vector field V on M, let n ) denote the local 

coordinates of a point p G M and (Vo(x), Vi(x), . . . , V n (x)) the covariant components of V at 
the same point. Thus, in case for any point p G M one has Vo(x) 7^ and V^(x) = dS{x)/dx tl 
for some S G C r (U p , R), r > 2 and U p denoting the local chart, we can apply Theorem 1' 
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(or the analogous for a metric having a different signature) to get the desired coordinate 
transformation on the generic chart U p . 

A local version of this global result is commonly quoted in General Relativity textbooks. 



IV. STOCHASTIC QUANTIZATION OF THE FREE RELATIVISTIC 

SPINLESS PARTICLE 

Nelson's stochastic kinematics for a spinless quantum particle is obtained by promoting 
the classical configuration variables to a Markov diffusion with constant coefficient equal to 
h/m, h being reduced Planck's constant and m the mass of the particle: 

dq(t) =b(q(t),t)dt+ (-Y dw(t) , *e[0,+oo). (4.1) 

\mj 

As briefly exposed in Sec. I, the hydrodynamical (or Eulerian) picture of the diffusion 
process is given, under appropriate regularity conditions, in terms of the time dependent 
density p and the current velocity v. The dynamics is then assigned by imposing that the 
motion extremizes the mean (regularized) classical action written as a functional of both 
fields |[L0|| , PH|| . In the case of a particle subjected to a scalar potential V, the time evolution 
of p and v is described by two coupled partial differential equations: 

d tP + V-( P v)=0, ^ 

in which v must be irrotational in all points where p is positive. Introducing therefore a 
scalar field S such that v = VS/m and performing the standard substitution 

^^«p{is ( ,,«)}, 
one can transform (|4.2| ) into the Schrodinger wave-equation 



dip r h 2 



— V 2 + mV 

2m 



iP . (4.3) 



Our objective is now to extend such a procedure to the quantization of a free relativistic 
spinless particle. This will be achieved by exploiting a suitable comoving coordinate system 
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and starting, as previously done ||12|| , from the classical Einstein action. We conjecture that 
the quantum evolution equation can be written also in the relativistic situation in terms of a 
couple of fields (p, v); they are expected to satisfy a covariant hydrodynamic equation which 
reduces to ( |4.2| ) (for V = 0) in the limit \\v\\ -C c. The relevant co moving coordinate system 
will be that induced by the four-velocity field V associated to v. Notice that the field V 
enters at this level as an unknown; the fact that it fulfills both assumptions of Theorem 1' 
will be verified a posteriori. To assign the kinematics we also need to introduce an evolution 
parameter independent of the particular choice of the coordinate A on 70'. To this end let 
as before /(A) denote the length of the arc (O', A) of 70' and consider a temporal parameter 
t defined by t := /(A)/c. By construction we have 



dl(\f = d£°d£ = -/(A) d\ 2 = -c 2 dt 2 , dt = v ~ dl ^ = - \/-£oo(£°) d£° . (4.4) 

c c v 

The physical meaning of this choice is that t coincides with the proper time of the comoving 
observer having 70' as its world-line. In the following we will also identify O = O'. 

We formulate the stochastic kinematics of a free relativistic spinless particle as follows: 

hi) there exists a timelike four- velocity field V on M 4 , satisfying the assumptions of 

Theorem 1', such that the evolution of the position of the particle with respect 

to the parameter t is described by a Markov diffusion in the comoving coordinate 

system associated to V; the diffusion coefficient is equal to H/m, m being the rest 

mass of the particle; 
/12) the current velocity of such a diffusion is equal to zero. 

Denoting by q l , i = 1,2,3, the spatial comoving coordinates of the particle on E , the 
time development of the Markov process is fully specified by writing 

dq\t) = Pi(q(t),t)dt+ f-Y G l k (q(t))dw k (t) , (4.5) 

where w is a standard Wiener process on R 3 and G\.G k i = g % \ g 1 ^ = o % i being the components 
of the (positive-definite) metric tensor of So (Corollary 2'). We will also assume that 
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the functions f3\ are sufficiently smooth to allow the existence of the so-called "backward 
representation" of the process pj||LT|| , namely 

dq\t) = (3l(q(t),t)dt+ ^ Gi(q(t))dv£(t) , (4.6) 

where is a reversed standard Wiener process. 

By introducing the infinitesimal forward and "backward" increments 



d+q^t) := q^t + dt) - q^t) 
d'q^t) := q\t) -q\t-dt) , 

we can construct the two four-vectors 

d+f := {d\, d+q^t)) 
dTp := (dX, d~q\t) ) 

in terms of which we define the mean regularized invariant 



(4.7) 



(4.8) 



E{d + Cd~^} = E{d + q\t) d- qi (t)} - c 2 dt 2 

= Eiv^qit^v^^qit^dt 2 - c 2 dt 2 . (4.9) 

Here eq. ( |4.4| ) has been used and v+, V- are to so-called "invariant drifts" introduced in 



Refs. [Kj and JT6|; they are vector fields on Eo defined by 



i ._ at h_ Q jkYi 

T* fc being the Christoffel symbols associated to gij (metric connection). 

As usual in Stochastic Mechanics, we introduce the current velocity (3 and the osmotic 
velocity u by combining the drifts: 

P ■= §(«+ + «-) = W+ + P-) (410) 



u 



\{v + -vJ) 



2 



It is convenient to introduce a smooth invariant measure on the manifold So and define the 
relative covariant probability density of the process, p(q,t), through the equality 
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E{F(q(t))} := J^F(q) p(q,t) J\g~\d 3 q , \g\ := det[^] , 

where F(q) is any Lebesgue integrable function of the process. It is well known [RJ that, in 
case p is strictly positive, (3, u and p are connected by the following relations: 

u* = — VMnp (4.11a) 
2m 

dtp + Viip?) = 0, (4.11b) 

where V« denotes the covariant derivative along direction q± (recall that by construction 
^%9jk — Vi, j, k). For the extension to the case when p has zeroes see [llj. According to 
/i2), we should also put (3 equal to zero (so that (3± become time independent); we prefer to do 
this later, in order to show that this choice is naturally enforced by consistency requirements. 

Thus, in analogy with the classical case, we consider, for arbitrary events a and h in M 4 , 
the following stochastic version of Einstein's action: 

A M := -mc J h ^E{-d+^d~^} = -mc 2 j*" Jl-E^^^dt , (4.12) 



which can be written as a functional of (3 and p exploiting (|4.10[) and ( |4. 1 la| ) . We get 



A [aM = - mc 2 f b dt (l - \ [ {f3 2 -u 2 )p^\ d 3 q) * , (4. 13) 

J t a \ C J / 

with the convention (3 2 := f3 l f3i = gij(3 l f3^ and similarly for u 2 . We recall here that, since the 



configurational manifold has zero intrinsic curvature, no Pauli-DeWitt curvature term fl7 
has to be considered in the expression of the stochastic action functional. 

Following a variational strategy which is analogous to that exploited in non relativistic 



Stochastic Mechanics ||15|| , we can now extremize the action with respect to independent 



variations Sf3 and dp. Since p and j3 are linked by the continuity equation ( |4.11b| ), this may 



be done by introducing a Lagrangian multiplier x, so that we construct the modified action 
A[a,b] ■= A [aM + f 6 dt f X [dtP + Vi{p(3 1 )} J\g~\d 3 q . (4.14) 

Jt„ JR 3 



/R3 

We require that 

14 



6vA[ aib ] = o(6p) 

8pA[ a ,b] = o(8p) 

for independent variations 5(3 and 5p having compact support on S x [t ,t&]. A straigth- 
forward calculation yields the two variational equations 

m(3 k 



'1-E 



V fc x 



k — X . 2 . 3 ; 



i 

2 



m f3 2 



+ 



m 



i 



2 2m 



We eliminate x by taking the gradient of the second equation; we are left with the following 
three-dimensional dynamical equation: 



dt 



'l-E 



n 



(4.15) 



2m 



Uk 



being VjV 1 the ordinary three-dimensional Laplace-Beltrami operator. 

In order to get a covariant expression of ( f4.15|) , we need a fourth equation which in the 
classical case is obtained by imposing that the action be stationary with respect to variations 
of the evolution parameter and which represents the conservation of energy. In this way 
the Hamiltonian is constructed as the conjugate momentum with respect to time. Let us 
then consider t as a new dynamical variable t = t(r), function of some invariant auxiliary 
parameter r. Since the Lagrangian does not depend explicitly on t, the requirement of 
stationarity of ( |4.12| ) with respect to smooth variations of t having compact support in 
[r a ,T b ], that is 

<M[*(*-.),tfo)] = °( St ) ' 
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immediately yields 



d 
If 



mc 



. 



Il-E 



(4.16) 



Alternatively we can also identify the Hamiltonian with a first integral of the motion. If 
(p*, denotes a given solution of ( |4.15| ) and ( [4.11b| ), then in the total variation of A only 
the boundary term survive: 



> --» ■..-./, / dt^ I X$P*\/\g\ d 3 q ■ 



8 nj . rA\, 



dt Jr 3 

We particularize now the variations in the following way: 

Sp* = p*St + o(5t) 
Sp, = $Jt + o(5t) , 



(4.17) 



so that expression ( |4.17|) becomes 



t„ Ot JR 3 V 



(4.18) 



On the other hand, being the dynamics time independent one can also write 

+ o(St), (4.19) 
so that, comparing ( 4.18 ) and ( 4.19| ) and taking the limit for St going to zero, we find 

(4.20) 



d 
It 







for any solution (p, (3) of ( |4. 15| ) and (fOTFj) (we note that L = L in this case). From this, 
continuity equation and an integration by parts give 

d mc 2 + E{mu 2 } 



dt 







(4.21) 



As a consequence, we can verify that a natural assumption in order the energy be properly 
defined is that the current velocity (3 in the comoving coordinates be equal to zero. 
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Under this assumption, the density of the process is made time independent so that 
dE{u 2 }/dt = by ( |4.11a| ) and ( |4 . 2 l| j becomes identical to ( |4.16| ). This discussion shows 
that hypothesis /12) on the kinematics is in fact crucial in formulating a mathematically and 
physically consistent theory. 

Putting explicitely hereafter (3 = and defining the constant functional 

1 

,2 
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1 + 



lip] 



(4.22) 



we get the following equations of motion in the comoving coordinates: 



M V^ Uk + (u&^Uk = 



f^mc 2 = 0. 



(4.23) 



It is apparent that for every fixed solution u of ( |4.23| ) the quantity 7 acts as a strictly positive 
constant; in fact, it drops out from the equations and we can solve the equivalent problem 
given by 



2 m 



(4.24) 



mc 







where the original comoving coordinates £ l ) are reintroduced (recall that Vj = V 
V ?i , i = 1,2,3). 



In order to express ( |4.24| ) in a covariant form, one has to introduce objects with the 

correct transformation properties. First of all, let us define a real field p by 

\ -1 



P(0 ■= N p(q) , 



Af 



-A 



-^oo(e )^ 



(4.25) 



where the constant Af has been chosen by demanding that p(£) is a normalized probability 
density on a fixed four-dimensional domain Q := [—A, A] x R 3 , A e R + : 

/ P(0 y/^ft = [ A M -9oo(Z°) d£° I J\g~\d 3 qp(q) = 1 , (4.26) 

with \g^ u \ := detfg^]. By construction p is a density with respect to the invariant measure 
on fi, so that it transforms scalar. 
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As a second step, we introduce also the two four-dimensional fields with components 

V» = I i ° , = A , = (O.t**) • (4-27) 

-#Qo(£°) 



Both [/ and V are four-vectors, due to different reasons. For the latter the four-vector 
character directly follows from the tensorial relation: 

U» = — VHnp , (4.28) 
2m 



which we obtain from ( |4.11a| ) and ( [4.25|) . The former is instead the four- velocity field 



corresponding to the current velocity vector (3 written in comoving coordinates (where the 
spatial part j3 is identically vanishing). It can be noticed that V M satisfy itself a relation 
similar to flPgQ : 

=, I V "S , S(0 := -mc [ e df \/-goo(?) ■ (4.29) 

Tfl J — oo 

It is important to recognize that eq. ( |4.29| ) automatically defines 5 as a scalar field. In 
all generality, the metric function goo(£°) must be negative definite, normalized according 
to fl4.26|) and sufficiently smooth to guarantee the existence of the primitive in ( 4.29 ). The 



simplest choice is obviously to take <?oo(£°) = — 1, or, physically, to adopt the proper time 
of the comoving observer along 70 as a coordinate, g?£° := cdt. In this case, the function S 
simply reduces to 

5(0 = -mc£° , e = ct. (4.30) 

Maintaining the above made positions and exploiting the fact that, by ( |4.27| ), U° = 0, 
V" 1 = 0, i = 1, 2, 3, we can now easily rewrite eqs. ( |4.24j ) as follows 

-(. 



A 2 

\2m 



. (4.31) 



The continuity equation can also be rewritten, in terms of a four- dimensional divergence: 

V M [pV-"] = 0. (4.32) 
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Now, by means of the same change of variables exploited to get the Schrodinger equation 
( |4.3p from its hydrodynamical counterpart (|4.2|) , it can be easily realized that ( |4.31| ) and 
( |4.32| ) correspond respectively to the real and imaginary part of a stationary four- dimensional 
Schrodinger equation, namely 

-|^ V <W = ( 4 -33) 

where 

fa) ■= VP(0exp{^(0}. (4.34) 

We remark at this point that, had we maintained 7 (eq. ( |4.22|) ) in the derivation, then 
we would have found an apparently non linear equation, which is quoted in |TJ]. We observe 
also that, as usual in Stochastic Mechanics, the eigenvalue [i is equal to the mean value of 
an energy. In the particular case one can show by direct computation that 

fi= I \\mV^ + \mU^\p^)^\d^ = -\mc 2 + \ m E{u 2 } . (4.35) 

We notice that, since u depends on the quantum state through p, the Schrodinger-like 



equation (|4.33|) exhibits a spectrum of energy levels. On the other hand, one can also observe 



that the ratio of the two contributions in ( |4.35|) is 
E{u 2 } h 2 



J 3 V^lnpV M lnp^jd 3 g , (4.36) 



c 2 4 m?c 2 

so that, avoiding extremely sharp densities, which would be patological in the free case, one 
can put with good approximation 

E{u 2 }/c 2 < 1 . (4.37) 

This approximation acquires a clear physical significance if the right hand side of ( [4.36| ) is 
interpreted as the (squared) ratio between the Compton wavelength h/mc and the typical 
length scale associated to the density p. Converting to energy units, one finds that the 
typical energy involved in the description has to be much less than 2mc 2 . This is precisely 
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the energy threshold under which a relativistic single particle theory is expected to hold. 
Turning to the four-dimensional Schrodinger equation ( f4.33| ), validity of ( |4.37| ) means that 
energy levels become degenerate and consequently the spectrum collapses to 

1 n 



mc 



Incidentally, this relation, also called " Feynman ansatz" , represents a widely accepted work- 
ing assumption in the literature |J, Equation ( |4.33| ) reduces then to the Klein-Gordon 
equation in comoving coordinates: 



m 2 c 2 



Tr 



4>(0 = o 



ten 



(4.38) 



Transforming to cartesian coordinates {x^} on M 4 , we find from this the familiar expression 

d d 



m 2 c 2 



dx^ dx u 



Tr 



[x 



o , x e®- 1 (ft) , 



(4.39) 



where, from ( 4.34|) , 



0(x) := V5W«P{^W} = >»«P{^«««))}- 

Since the parameter A which fixes the domain fl is arbitrary, the region where the equation 
is defined can be sufficiently large to contain any arbitrarily chosen rectangle in M 4 . 



V. STOCHASTIC INTERPRETATION OF DIFFERENT SUBSETS OF 
SOLUTIONS TO KLEIN-GORDON EQUATION AND 
NON RELATIVISTIC LIMIT 

The quantization procedure which we have just outlined naturally selects a particular 
subset of the solutions to Klein-Gordon equation. In fact we observe that, in comoving 
coordinates, the conserved four-current density at an arbitrary point x G M 4 , defined as 

>(0 := mp(0 V"(t) = p(0V^(0 , £ = Hx) , (5.1) 

has the following particular structure: 
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J°(0= m ° m \, , ^(0 = 0- (5.2) 
-9oo (6°) 



We also notice, that the conservation law associated to (|5.2|), namely 



= 



is nothing but a rewriting of the continuity equation (|4.32|) for the probability density of the 
stationary diffusion q(t). 

For a generic solution of the Klein-Gordon equation in the inertial frame let us corre- 
spondly introduce the conserved four-current density by 

J„{x) := |0(x)| 2 ^S(x) = YiifW-tW) ■ ( 5 ' 3 ) 

Denoting by V(x) the four- vector with covariant components V^{x) = d ll S{x)/m, we have 
by construction 

W(0 = ^(x)V»(x) } t = *{x). (5.4) 

We then choose in a neighborhood of £ normal coordinates := (£°, £l, As observed 

in Sec. Ill, we get 

= A-'AxW^x) , (5.5) 

where Vf denote the normal components of V and is, up to a spatial rotation, the Lorentz 
boost associated to the three-dimensional velocity 

u*:=c^, 1 = 1,2,3. (5.6) 

Thus we can conclude that the comoving coordinate transformation $, represented in 
normal coordinates in the neighborhood of a point x by $ = $*, acts on four- vectors as 
a Lorentz boost associated to the three- velocity v(x) (which is uniquely determined by the 
four- velocity V(x)). 

We also observe that since J M is itself a four-vector, then 
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= as(£K°(0, (5.7) 

with (from Q) 

J°(0 = mcp(0 , (5.8) 

Recalling the transformation properties of p and expliciting the Lorentz boost, one finds 
from 

_ mcp t _ rmfp 

so that 



\/x e . (5.10) 



J° > 

Thus we can immediately recognize that a necessary and sufficient condition for a solution 
to Klein-Gordon equation corresponds to a one-particle diffusion satisfying the variational 
principle in the comoving coordinates is that conditions fl5.IGp hold. 

In fact, necessity was just proven. Conversely, given a normalized solution of the Klein- 
Gordon equation with associated four-current density J, sufficiency immediately comes by 
putting, in all points where J° ^ 0, 

v % :=c 3 — , p:=0*0. (5.11) 



We have from this 



and consequently 



J» = -p T = Aj* , J- = ^JL = A;j? , (5.12) 



which proves the assertion. 
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Some remarks are now in order. First of all, we notice that the set of solutions restricted 



solutions (in a suitable limit of non normalized 0). Note also that the positivity constraint 
on J° consistently allows to interpret it as a conserved density in hydrodynamical sense. 

Furthermore, this stochastic quantization procedure can be extended in an obvious way 
to the description of a beam of identical (spinless) non interacting particles: in this case, 
one would simply recover positive energy solutions fulfilling both ( |5.10|) but such that 



N denoting the total number of particles. (This is trivially obtained by substituting the non 
normalized density pN in place of p at the beginning of our derivation). 

A less obvious fact is that it is possible to give a fully probabilistic interpretation also 
to solutions with negative energy. The idea can be sketched as follows. Let us start from 
a R 3 -valued diffusion having zero current velocity and satisfying the stationary stochastic 
differential equation 



dq(t) = /3+(q(t))dt + v*dw(t) = u(q(t))dt + u*dw(t) , t G [0, T] , dt > , (5.13) 



with u denoting, as usual, the osmotic velocity and T > arbitrary. Let us now think of 
a diffusion with a "specular" time evolution. In order to have an example of a (classical) 
phenomenon that would be described by means of such a "specular" diffusion, one could 
imagine a cloud of brownian particles which, starting from an initial spreaded spatial dis- 
tribution, would concentrate after some interval of time in a small region. This is a "rare" 
event to which usual probabilistic models would give zero probability, but which is not, of 
course, in principle impossible. To get a mathematical description in terms of "forward" 
differentials, let us firstly introduce time-reversal in the usual way: 



by conditions ( |5.10| ) is not empty since it contains at least all positive frequency plane-wave 





(5.14) 
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where 5 is a fixed positive time increment. We have immediately from (|5.14 ) 



dq\t' - 5) := q\t') - q'(t' - 5) = - dq(t) =, - (q(t + 5)- q(t)) . (5.15) 

But dq(t) is assigned from eq. ( pM3Q ; in particular, we can use the corresponding so-called 
"backward representation" to write 

dq'(t'-5) = -dq(t)= -{-u(q(jt))dt + v* dw*(t)} , (5.16) 

where we have exploited the relation (3- = — f3 + holding in this particular stationary case 
flQfl , and where is a standard reversed Wiener process, that is a process which has all 
properties of an usual Wiener process except from the fact that its increments are indepen- 
dent of the cr-algebra generated by the future of w*{t). Let us introduce 

w'{t') - w'{t' - 5) := w*{t) -w*(t + 5) . (5.17) 

From the definitions and the properties of w*, one can immediately see that w'(t') is a 
standard Wiener process (which, loosely speaking, "goes forward in time"). By inserting 
eq. ( |5.17| ) in eq. ( |5.16| ), we are able to write, after standard manipulations, a stochastic 



differential equation for the diffusion q', namely 

dq'(-t) =u(q(-t))dt + v5dw'(-t) , -t E [-T, 0] , dt > . (5.18) 

Eq. ( |5.18|) provides a possible forward description of what we have called a "specular dif- 
fusion". Now, all of these considerations can be extended to the stationary kinematics in 
comoving coordinates. We can thus conclude that to every solution of the Klein-Gordon 
equation with positive J° there corresponds a solution where t is replaced by —t. The first 
one is stochastically represented by a stationary diffusion in the comoving coordinates (with 
respect to the invariant parameter t), while the second one corresponds to the "rare event" 
of a "specular diffusion". Changing the parameter t into —t simply leads to a change of sign 
in S, as one can easily see recalling eq. ( |4.30| ) and the fact that 5* is a scalar by construction. 
Then J° becomes negative and consequently 
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is also negative at every point. This furnishes us with the desired interpretation. 



Having clarified how different subsets of solutions to Klein-Gordon equation can receive 
different stochastic interpretation, we are now in position to study the non relativistic limit 
described by Schrodinger equation. Since the latter represents the evolution equation for the 
quantum state of a single standard spinless particle, we need to consider only those solutions 
which satisfy conditions ( |5.1C| ); in addition, we must also select those for which ||u|| <C c. 
We expect that these last solutions approximately solve dynamical equations equivalent to 
the Schrodinger one. 

The simplest starting point is the Klein-Gordon equation written in cartesian coordinates 
d d 



m 2 c 2 



4>{x) 







(j){x) 



dx^ dx u h' 

being p(x) a positive normalized probability density over M 4 . If we define 

1 d 



p(x) exp j— *S r (rr) ^- , (5.19) 



m dx^ 



S(x) 



(5.20) 



then, in agreement also with ( 5.11 ), there exists a velocity field v on R 3 such that 



1 _ vi 



1 _ VI 



(5.21) 



We also construct the four- dimensional osmotic velocity by 

h d 



u„ 



2m dx^ 



\np(x) . 



(5.22) 



The Klein-Gordon equation (|5.19|) can be now expressed in terms of U and V: by following, in 
reversed order, the same steps leading from eqs. (|4.31|) and (|4.32|) to eq. (|4.33| ), the real part 
and the gradient of the imaginary part constitute a pair of four-dimensional hydrodynamical 
equations: 
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Jl 

dx 



\2m 'I dx^ dx" dx^ J dx^ 

'-(pv 11 ) = o. 



(5.23) 



Starting from ( 5.23| ), we now claim that the non relativistic limit is recovered, after 
substituting V° and V % in terms of ( 5.21 ), by taking the limit \\v\\/c <C 1 in the resulting 
expression, having care to treat as negligible all terms which are multiplied by h 2 /c and 

h 2 /c 2 . 

The proof goes as follows: let us denote by u the spatial part of U and by v the three- 
dimensional velocity field associated to V, and let t be the time parameter in the inertial 
frame associated with {x^ 1 }. We put, as usual, dx° := cdt. Letting ||u||/c <C 1, the four- 
dimensional continuity equation can immediately be rewritten as 

dtp + V • (pv) = . 

The second hydrodynamical equation is now considered in two steps. Separation of the 
spatial part leads to the following three-dimensional equation: 



V 2 u 

2m 



h d 2 u 
2^~ch^ 2 



dx - 



+ 



■ V 



i rr_ 



i v 



+ 







1 _ dx° 



1 V 

r _q 



. (5.24) 



If, as before, we rewrite dx° = cdt and take into account the explicit expression of U fl5.22| ), 
we can neglect the second and fourth term, both of which depend on h 2 /c 2 . In the limit 
||u|| <C c, Eq. ( |5.24 ) becomes then 



h 

— V 2 m + (u • V) u 
2m 



/ ^ . \ dv 

+ ( v .y) v + 

at 



. 



This, together with the continuity equation, is nothing but the hydrodinamical form of 



Schrodinger equation, which we have also quoted in Q4.2 ). The temporal component of the 
dynamical equation ( |5.23j ) gives rise instead to the following one- dimensional relation: 



n 

2m 



—V 2 U n 



+ 



1 - % 

cr 



1 - ^ 







■£_ dx° 



1 - ^ 



. (5.25) 
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The last two terms are readily recognized to be neglegible in the limit \\v \\/c <C 1. Recalling 
the explicit form of Uq coming from ( |5.22| ), we see that the other terms depend on h 2 /c. 



This proves the assertion. 

It is worth to note that the quantum non relativistic limit appears as an approximation 
of the correct relativistic dynamical equations from a twofold point of view: firstly, the 
approximation \\v\\ C c corresponding to the classical non relativistic limit is done; secondly, 
terms depending on h 2 /c are also neglected. This fact is perhaps a little unexpected and 
conceptually intriguing. 

Furthermore, an important point of this approach needs to be stressed here, namely 
the fact that in the relativistic framework no diffusion process is supposed to exist for the 
inertial observer. A simple probabilistic interpretation is consequently lost in a generic 
inertial frame. This difficulty may be considered, in a sense, as the stochastic counterpart 
of non trivial questions which, in the framework of canonical quantization, are related to a 



consistent definition of observable quantities, such as position operator (see |7j and [19] for 
any detail). In our procedure, the reason for this fact is that, as observed above, the Wiener 
process has sample paths which can "go outside" the light cone, so that there exist inertial 
frames where such spacelike trajectories would appear as going forward and backward in 
time. 

On the other hand, such a peculiarity disappears in the non relativistic limit. As a con- 
sequence, the usual Nelson's quantization becomes consistent with the proposed procedure 
in a non relativistic setting. 

VI. DISCUSSION AND OUTLOOK 

The first observation to be done is that stochastic quantization of the free spinless rel- 
ativistic particle based on comoving (local or global) coordinates is very clean both from a 
mathematical and a physical point of view. In fact, on one side all probabilistic objects are 
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well defined; on the other side, once given stochastic kinematics in the comoving coordinates, 
dynamical equations naturally emerge from a stochastic version of the Einstein action, by 
means of the stationary action principle. No additional ad hoc assumptions are necessary, 
at variance with the path integral approach |0| and previous attempts within stochastic 
frameworks which we have already quoted. 

We also observe that the stochastic quantization procedure proposed in this work gives 
only to the subset of solutions to Klein-Gordon equation with positive energy and normalized 
four-current the physical meaning of representing the quantum evolution of one single free 
relativistic spinless particle. As is well known, this is a consistency property which is not 
easily achieved within standard canonical (first) quantization ||19|| . 

On the other side it may be interesting that one can give a probabilistic interpretation 
also to other solutions of Klein-Gordon equation, in particular to those with negative energy 
in terms of specular diffusions as " rare events" . 

It may also be worthwhile to insist on the fact that the comoving coordinates approach 
allows to more properly understand how to handle covariant and non covariant quanti- 
ties. Dynamical equations (which do not contain stochastic terms) are Lorentz covariant, 
while Markovian kinematics is not: in fact, the latter is covariant only with respect to 
reparametrizations on Eo and 70 

We point out that global comoving coordinates approach introduces a description of the 
motion of the particle which is well defined on the whole M 4 (this in particular provides the 
technical advantage of eliminating all spurious boundary terms which come out in the local 
derivation |f20|| ). Moreover, it directly leads to a dynamical theory which is covariant with 
respect to arbitrary changes of space-time coordinates. 

Indeed the results presented in this paper look suitable for the extension to the description 
of a quantum particle subjected to a gravitational (classical) field. In this case one can 
furthermore conceive that, in presence of strongly peaked gravitational fields, the degeneracy 
among the energy levels of the four- dimensional dynamical equation ( f4.33p breaks down. 



This is an aspect of the outlined stochastic quantization procedure which deserves further 
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investigation pi 



From another point of view, it would be also appealing both on mathematical and physi- 
cal grounds to consider a sort of second stochastic quantization, by reinterpreting the proba- 
bility density in the comoving coordinates as a physical scalar field in a space-time dependent 
random medium. This could in fact be done by exploiting suitable probabilistic techniques 
(see Ref. 231 for an excellent review of this subject). 
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